RKKY interactions and anomalous Hall effect in metallic rare-earth pyrochlores 
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Motivated by experiments on Pr2ir207, we consider metallic pyrochlore systems A2B2O7, where 
the A-sites are occupied by rare-earth local moments and the B-sites host 5d transition metal ions 
with itinerant strongly spin-orbit coupled electrons. Assuming non-Kramers doublets on the A- 
site, we derive the RKKY interaction between them mediated by the B-site itinerant electrons and 
find extended non-Heisenberg interactions. Analyzing a simplified model of the RKKY interaction, 
we uncover a local moment phase with coexisting spiral Ising-like magnetic dipolar and XY-like 
quadrupolar ordering. This state breaks time-reversal and lattice symmetries, and reconstructs the 
B-site electronic band structure, producing a Weyl Metallic phase with an intrinsic anomalous Hall 
effect and an undetectably small magnetization. We discuss implications of our results for Pr2lr207- 



The metallic pyrochlores, A2B2O7, with a rare-earth 
A-site ion and a 5d transition metal B-site ion, lie at the 
intersection of exciting recent developments in condensed 
matter physics. The rare-earth moments on the A-site 
pyrochlore sublattice could lead to quantum spin ice 
physics [H [2] from Ising anisotropy and geometric frus- 
tration, a feature they share with the well-studied insu- 
lating pyrochlore oxides. The strongly spin-orbit coupled 
5d conduction electrons, on the other hand, contain the 
seeds of topological phases like Weyl semimetals or topo- 
logical insulators [3HH] • The interplay of these two effects 
could pave the way for new emergent phenomena. The 
pyrochlore iridate Pr 2 Ir 2 07 provides an example of such 
a metallic frustrated system, with a significant Curie- 
Weiss temperature 9cw ~ — 20K, much larger than in 
the insulating compounds Pr 2 Sn 2 07 or Pr 2 Ti 2 07 [5HT2"]. 
For T < 1.7K, even when an applied magnetic field along 
the (111) direction is switched off, it exhibits a significant 
anomalous Hall effect (AHE) which grows upon cooling 
[T5HT5] . although, in contrast to Nd 2 Mo 2 07, the uni- 
form magnetization is undetectably small over a range 
of temperatures. [16] This raises the key issue of mecha- 
nisms underlying the AHE [17] in materials with geomet- 
ric frustration and strong spin-orbit coupling. 

For metallic pyrochlores, the RKKY interaction be- 
tween the local /-moments induced by the conduction 
electrons is expected to be important. When the /- 
moment is a non-Kramers doublet, such as in Pr 3+ , the 
pseudospin component t z along the local (111) direction 
carries a magnetic dipole moment and couples to the con- 
duction electron spin density, while the r x ' y components 
carry a quadrupole moment and couple to the electronic 
charge density. We show that the resulting RKKY cou- 
pling has a highly non-Heisenberg form and extends be- 
yond the nearest neighbor term. We propose and study 
a simplified such extended XXZ model for non-Kramers 
doublets using a variational analysis, finding a rich phase 
diagram, shown in Fig. [I] which includes incommensu- 
rate spiral 'coexistence' states (CE l5 CE 2 ), with wavevec- 
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FIG. 1: Variational phase diagram of the pseudospin Hamilto- 
nian in Eq.|2| relevant to RKKY coupled non-Kramers dou- 
blets as a function of anisotropy w and 8 — tan~ 1 (J2/Ji). 
w > 0: We find ferro-XY (FM-XY), antiferro-XY (AF-XY), 
where spins order in the plane transverse to the local (111) 
axes. These are time-reversal invariant states with modulated 
quadrupolar order, w < 0: We find Ising all-in/out (AIO), 
Ising two-in/out (2IO) and states with coexisting modulated 
Ising- XY orders (CEi,CE2). The CE2 state breaks time- 
reversal and lattice symmetries, producing an anomalous Hall 
effect for the spin-orbit coupled conduction electrons. 



tors (0, q, 7r) or (0, 0, q), having spatially modulated, mag- 
netic (dipolar) and quadrupolar order, i.e., 'magneto- 
quadrupolar supersolids'. We confirm that such a CE 2 
state also exists within a classical Monte Carlo simula- 
tion. Although the CE 2 state has no net local moment 
magnetization, it reconstructs the band structure of the 
spin-orbit coupled 5d-electrons, producing Weyl points 
as well as small Fermi pockets, leading to a measurable 
AHE and an extremely small conduction electron mag- 
netization [21] • We discuss possible implications and pre- 
dictions for Pr 2 Ir 2 07, for which our proposal appears to 
be distinct from the more widely discussed spin chirality 
scenario [TS1 [TSHIJ] . 
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Basic microscopies and RKKY coupling. — Since we 
are motivated by experiments on Pr 2 Ir 2 07, we briefly 
review the relevant microscopies of the local moments 
and conduction electrons in this material. The Pr 3+ 
ion is in a 4/ 2 configuration, with electron-electron in- 
teractions and strong spin-orbit coupling favoring a to- 
tal angular momentum state with J = 4. This nine- 
fold degenerate manifold is split by crystal fields with 
D%d point group symmetry around the Pr 3+ ion, aris- 
ing from a cage of eight O 2- ions stretched along the 
local (111) axis. The splitting is captured by an effec- 
tive time-reversal invariant crystal field Hamiltonian ex- 
pressed in terms of J=4 angular momentum operators, 
H ccl = ~aJ*+/3[J z (Jl+,P)+h .c .]+ 7 (j6+j6 ) ) with a>0 
and /9,7 < a, leading to a ground state non-Kramers 
doublet with a dominant \J Z — ±A) component. [2"T] Pro- 
jection to this low energy doublet allows us to define 
pseudospin-1/2 operators t with t z cx J z while t oc 
{J±,J Z }. Under time-reversal, t z — > — r 2 , transforming 
as the magnetic dipole moment, while r x,y are left invari- 
ant, transforming like a quadupole moment. 

For Ir 4+ electrons, strong spin-orbit coupling in the 
t2 g manifold leads to a half- filled effective j=l/2 band. 
The metallic character of Pr2lr207 suggests that a tight 
binding model which ignores strong electronic correla- 
tions would be an adequate starting point, 
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where c\ a {c^) denotes the electron creation (annihila- 
tion) operator at site i with the Kramers pseudospin in- 
dex a corresponding to j z =±1/2, and er = (a x ,a v ,a z ) 
are the Pauli matrices. We assume tij = t (i') for nearest 
(next-nearest) neighbors, and Vij ^ only for nearest- 
neighbors and is constrained by lattice symmetries. |22j 
How does the A-site non-Kramers doublet couple to 
the B-site conduction electrons? Unlike the usual Kondo 
coupling to a magnetic Kramers doublet, time-reversal 
invariance dictates that the Ising component of the A-site 
pseudospin r? at site j (which points along the local (111) 

axis) couples to spin density jf = <4a~f^ c i/3 °f electrons 
on the six neighboring B-sites, while the planar compo- 
nents of the pseudospin rp couple to the charge density 

n i — c ia C ia on the neighboring B-sites. Keeping j^Tj 
and TiiT^ terms, we find the symmetry allowed Kondo 
coupling Hab with three parameters ci, C2, C3. (See Sup- 
plementary Materials for details.) Integrating out the 
conduction electrons, we find that the resulting RKKY 
interaction has two important features which are insen- 
sitive to the details of the Ir band structure: (i) it allows 
for significant couplings beyond the nearest-neighbor in- 
teraction, but is negligible beyond the third neighbor; 
(ii) it is highly anisotropic in spin-space since t x,v and 
t z interact very differently with the spin-orbit coupled 
conduction band. A complete microscopic set of Kondo 



couplings should include rjc. 

i 7^ £, terms which we have omitted, and the full char- 
acterization may require multiple parameters. Rather 
than dealing with such a complex model, we would like 
to focus on general and robust features of the RKKY in- 
teractions. We therefore turn to a study of a simplified 
RKKY Hamiltonian which retains the two key features 
described above. 

Local Moment Model and Phase Diagram. — Moti- 
vated by our observation that the RKKY interaction be- 
tween local moments is highly anisotropic in spin space, 
and has beyond nearest-neighbor terms, we study a sim- 
plified model in the basis with a sublattice- dependent 
quantization axis along the local (111) direction, 



H = J r 



{l^w)Tl{r)r!,{r')+Tjr{r)-Tjr,{r') (2) 



Here J r r > = J\{J2) for nearest (next-nearest) neighbor 
sites, T z (r) is z-component of the pseudospin on sublat- 
tice s at site r, tj~ (r) denotes the transverse component 
of the pseudospin which lies in the local XY plane, and w 
quantifies the exchange anisotropy. We set J\ — J cos 9 
and J2 = J sin 9, and explore the phase diagram of this 
model as a function of (9,w). 

Treating the spins as classical unit vectors, we mini- 
mize the energy using a variational ansatz 



T s (r) = d s e 3 + yi^I TZe[(e 1 +te 2 )e^- r+ ^}, (3) 

where £1,2,3 form a triad of orthonormal vectors, so that 
|r s (r)| = 1. In the local coordinate system, this ansatz 
allows, (i) for d s — 0, a coplanar spiral with wavevector 
Q with spins in the (ei,e 2 ) plane, and, (ii) for d 2 s = 1, 
a collinear state with spins along ±£3. In the isotropic 
limit, w = 0, this ansatz recovers (0, 0, q) spirals [2"3"] . 
while in the Ising limit it allows for 2-in 2-out or all-in 
all-out states. The complete ground state phase diagram 
from this variational analysis is shown in Fig. [T] 

For w > 0, we find states where the spins lie in the local 
XY plane, forming phases like ferromagnetic XY, antifer- 
romagnetic XY, or degenerate XY versions of 2-in 2-out 
states. These XY ordered states for the non-Kramers ions 
do not break time-reversal symmetry but correspond to 
modulated quadrupolar orders, with no net quadrupole 
moment and hence may be termed 'antiferroquadrupolar' 
(AFQ) states. With increasing 9, the different colors in 
the AF-XY region represent states with distinct ordering 
wavevectors Q = (Qqn), (0,0,0), (00g). 

For w < 0, by contrast, we find in addition to well 
known states like the ferromagnetic Ising (all-in all-out), 
and Ising spin ice (2-in 2-out), large parameter regimes 
which support coexistence phases (CEi,CE2) with copla- 
nar order involving spatially modulated Ising (magnetic 
dipolar) and XY (quadrupolar) order. The CEi and CE 2 
states order at the wavevectors (0, q,ir) and (0,0, q) re- 
spectively (or their symmetry related momenta). Such 



CE states are 'magneto-quadrupolar' supersolids which 
break time-reversal symmetry (defined by t z — > — t z ), 
and most lattice symmetries. Below, we focus on the 
remarkable physical properties including the AHE of the 
CE2 state which is (i) robustly present in an unbiased nu- 
merical energy minimization using a simulated annealing 
approach [3D] , and (ii) stable to the addition of weak per- 
turbations such as (t^t^ + t~tJ) to the simple Hamil- 
tonian in Eq. [2j A more complete study including such 
terms which break the local U(l) symmetry will be dis- 
cussed elsewhere [23] . 

AHE and Magnetization in the CE2 state. — Once the 
CE2 state is stabilized on the A-site, it imprints effective 
spatially varying magnetic fields and chemical potentials 
on the B-site conduction electrons due to the Kondo- 
type coupling between them. For an ordering wave vector 
Q, this mixes electrons with wave vectors k and k + Q, 
leading to a Hamiltonian of the form 



Hg{k) = 



Ftb(fc) Hab(Q) 
Hab(-Q) H th {k + Q) 



(4) 



For the CE 2 state in our anisotropic J\-J-z pseudospin 
model, Q = (OOq), and time-reversal symmetry and all 
lattice symmetries except Ci z (W rotation along z direc- 
tion) are broken. Once these broken symmetries are in- 
herited by the resulting 5d electron band structure, one 
can argue for an intrinsic AHE. The Ci z symmetry leads 
to a yz = <j xz = since the current operator J z is invari- 
ant under C*2 Z rotation (x, y, z) — > (— x, —y, z) while the 
currents J y and J x change sign. However, a xy remains 
unchanged under this rotation and can thus be nonzero. 

To explicitly compute the Hall conductivity, we assume 
the tight-binding parameterization for the Ir 5d elec- 
trons in Eq.([T]), as relevant to P^h^Oy, choosing t = 1, 
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-0.1, and \vij\ = 0.2. These values are close to those 



determined from the Slater-Koster parameters. [22] Fig{2] 
shows an example of the reconstructed B-site band struc- 
ture based on Eq. Q . (Red line indicates the Fermi level. ) 
We consider the A-site incommensurate CE2 local mo- 
ment order with wave vector Q ss 1.27r(001), which is ap- 
propriate for 6 — 1.72-K, w = —0.5 in Eq.Q (* marked in 
Figjl]), and the Kondo coupling constant c/t = 0.1 where 
c 1 = 0, C2 = C3 = c in Hab(Q)- (See Supplementary ma- 
terial) The reconstructed band structure generates both 
Weyl points and Fermi pockets near r, V + Q. There are 
four pairs of Weyl points in total (two pairs near r, T + Q 
and other two pairs at their C'2 Z symmetry related points) 
and the inset of Fig |2] shows one of those Weyl points and 
Fermi pockets near k sa (0, 0.127T, 0.027r). 

Figj3] shows the explicitly calculated Hall conductiv- 
ity using the Kubo formula as a function of Kondo cou- 
pling constant c, based on the reconstructed B-site band 
structures. The strength of Kondo coupling c determines 
the magnitude of the AHE response. a xy initially in- 
creases with increasing c, acquiring contributions from 
small electron-like and hole-like Fermi pockets as well as 
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FIG. 2: Reconstructed 5d conduction electron band struc- 
ture (based on Eq.Q with the energy units of t) in the pres- 
ence of Kondo coupling constant c/t = 0.1 and the A-site 
incommensurate CE2 local moment order with wave vector 
Q ss 1.27r(001). The inset shows one Weyl point and Fermi 
pockets near fc « (0, 0.12-7T, 0.027r) (the horizontal axis repre- 
sents k along k — (k, k + 0.12-7T, 0.027r)). Its pair is located 
at fc ss (O.O671-, 0.18-7T, 0). Such a pair of Weyl points is near F 
point and there is another pair near T + Q. Furthermore, ad- 
ditional two pairs exist at their C2 Z symmetry related points. 
Overall, there are four pairs of Weyl points. 



four pairs of Weyl points induced by the CE 2 order. For 
c/t > 0.7, a band gap opens up and a xy = 0. 

On symmetry grounds, the AHE must be accompanied 
by a nonzero uniform magnetization, with M x — M y = 
but M z ^ 0. Remarkably, although the CE2 state has 
no net magnetization from the local moments, we find 
that it induces a nonzero magnetization M z for the 5<i 

electrons, where M M = ^QZi Y, a /3 CiaT^ *?)- For sma11 
c, the net magnetization M z gets larger, proportional to 
the density of states (DOS) near the Fermi level, with the 
same sign and a trend qualitatively similar to a xy . For 
large c/t > 0.7, we find that M z ^ although a xy = 0, 
signalling a magnetized band insulator. 

This dichotomy of a large a xy but a small M can 

be argued for as follows. On dimensional grounds, the 

2 

AHE signal <j xy ~ j^Ak, where the momentum scale Afc 
must be induced by an effective 'internal magnetic field' 
B lnt due to the spontaneous breaking of time-reversal 
symmetry 25, 26J. If the t x,v order in the CE2 state re- 
constructs the band structure to produce a small Fermi 
pocket, with a Fermi wavevector Uf and an effective mass 
to*, we expect Afc ~ B lnt m* /Uf resulting in a large a xy 
due to a small hp Fermi pocket, while the magnetization 
M ~ B lnt m*kF stays small. The Weyl point contribution 
from splitting a quadratic band touching, with a curva- 
ture m** , on the other hand leads to Ak — Vto** _B int [2"T]. 
and a magnetization M ~ B mt , so that a small B lnt again 
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FIG. 3: Hall conductivity a xy and magnetization M z of 5d 
conduction electrons as a function of Kondo coupling con- 
stant c based on the reconstructed band structure. (See Fig[2] 
for the reconstructed band structure at c/t = 0.1) For c = 0, 
the conduction band is in a semi-metal phase (SM) with a 
quadratic band touching at the F point. For < c/t < 0.7, 
the CE2 state reconstructs the bands generating both Fermi 
pockets and pairs of Weyl points, resulting in a Weyl-metallic 
phase (WM), with a significant nonzero a xy which has con- 
tributions from Weyl points as well as Fermi surfaces. For 
c/t > 0.7, the conduction electrons form a gapped normal 
insulating phase (NI). The magnetization M z has the same 
sign as <j xy over a wide range of couplings c, although it is 
extremely small in magnitude, and is nonzero even in the NI. 



leads to a large a xy and a small M. The presence of both 
contributions would lead to a non- linear <j xy (M). 

Application to Pr 2 /r 2 07: The CE2 state in our sim- 
ple model of RKKY coupled non-Kramers ions and spin- 
orbit coupled conduction electrons captures a key as- 
pect of the experimental data on Pr 2 Ir 2 07: an AHE 
accompanied by a negligible magnetization. However, 
since our ordering wavevector is along the (001) di- 
rection, the spontaneous AHE is produced in the xy- 
plane, whereas the spontaneous AHE in Pr 2 Ir 2 07 is seen 
for fields along the (111) direction. This discrepancy 
may be resolved if such coexisting spiral order had a 
wavevector along the (111) direction, or if our CE order 
gave way to a multimode spiral formed by superposing 
(0, 0, q), (0, q, 0), (q, 0, 0) spirals while preserving C 3 rota- 
tion along (111) direction but breaking all other lattice 
symmetries. The Hall conductivity in that case will nat- 
urally be in the plane perpendicular to (111), in agree- 
ment with experiment. We estimate the magnitude of the 
AHE (at c/t = 0.1, which is c « 5meV) in our CE2 state 
to be a xy re -^(Q^cm" 1 ) and M z - 0.002 (/x B /Ir). 
This is consistent with the magnitude of the experimen- 
tally measured Hall signal, and the absence of any mea- 
surable magnetization over a range of temperatures. At 
lower temperatures, experiments detect a nonzero rem- 
nant magnetization, which may arise in our model from 
the feedback of the magnetized conduction electrons on 
the local moments, an effect we have not taken into ac- 



count. 

The Curie- Weiss temperature is sensitive to J2/J1, 
changing sign as a function of 9 even within the CE 2 
state. For couplings c/t ~ 0.1, we can obtain 9cw to be 
of the right sign and magnitude, 8cw ~ —20K as ob- 
served experimentally, for parameter values closer to the 
all-in all-out phase boundary. 

The CE 2 state should exhibit two distinct thermal 
transitions, associated with the onset of Ising and XY 
orders, with this splitting being smaller for a weaker 
anisotropy w. However, this simple expectation gets con- 
founded by two issues, (i) Terms which we have omitted, 
such as (t^t^ + t//tJ), will break the (staggered) U(l) 
invariance of the XY terms, and modify the 3d-XY uni- 
versality class of the quadrupolar ordering transition, (ii) 
Oxygen vacancy defects will produce strong electric fields 
which break the point group symmetry around the 
p r 3+ j on ^ splitting the non-Kramers doublet due to extra 
random terms AH CC { cx {J z , J±} ~ in the effective 
crystal field Hamiltonian. This would lead to a time- 
reversal invariant strong random field on t x ' v . These 
effects might conspire smear or destroy the XY ordering 
transition. However, the time-reversal symmetry break- 
ing Ising transition is expected to survive, which would 
be consistent with the single specific heat "peak" seen at 
the onset of the spontaneous AHE in Pr 2 Ir 2 07. Monte 
Carlo studies of the thermal properties and disorder ef- 
fects will be discussed elsewhere 24] . 

The most direct evidence for the scenario advocated 
here would be a probe which can detect the modulated 
Ising order using neutron diffraction. Landau theory ar- 
guments predict that entering such a coexistence phase 
in a clean system would also lead to a weak charge den- 
sity wave of the Ir electrons at the ordering wavevector 
of the spiral, from a term cx T x ' v (q)p(—q), and which we 
find small but finite in our calculations; such charge or- 
der may be weak but can be probed, in principle, using 
X-ray diffraction. In the presence of an induced charge 
order, Landau theory arguments also predict a nonzero ri- 
band magnetization, from a term cx Mt z (q)p{—q), which 
is indeed present as discussed above. Quantum oscilla- 
tion measurements to detect the ordering-induced Fermi 
pockets are desirable, but might be difficult due to the 
significant AHE and the small size of the pockets. 

Conclusion. — In conclusion, we have proposed a mech- 
anism of intrinsic AHE in metallic pyrochlore systems, 
such as Pr 2 Ir 2 07, arising from spiral order of local mo- 
ments driven by their extended anisotropic RKKY ex- 
change interactions, and the resulting reconstruction of 
the electronic band structure to form small Fermi pock- 
ets and pairs of Weyl points. This ordering could occur 
proximate to an all-in-all-out state of the local moments, 
and appears to be distinct from previously proposed spin- 
chirality scenarios for the AHE in Pr 2 Ir 2 07. 
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Supplementary Material 
Local coordinates and primitive vectors in A2B2O7 

First, we define the local coordinates for every A-site 
pyrochlore lattice. The local spin quantization axis (z) is 
always pointing towards the center of each tetrahedron 
and x, y axes are defined on its basal plane. Table[I]shows 
their local coordinates frames depending on the four sub- 
lattices of the pyrochlore lattice. The A-site pseudospin 
in global coordinates, say f s (r), can be represented as 



,(r) 



(r)z s 



(5) 



where r s (r) = (t® ( r ) > T s ( r ) > T t ( r )j is the pseudospin at 
site r and sublattice s defined in the local coordinates. 
Both A and B-sites form pyrochlore lattice and each can 
be viewed as FCC lattice with a four-site basis. We de- 
fine the primitive vectors (h) and the position of four 
sublattices (Ai,Bi) for A and B-sites in Table [TT] 
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3 


Xi 


71 ( 0l1 ) 


75( 011 ) 


75( 011 ) 


75(° n ) 


Vi 


7f(2H) 




7f( 211 ) 


7f( 211 ) 


Si 


Tf^ 111 ) 


71 


7f( IlT ) 


7i( ni ) 



TABLE I: Local coordinate frames for the four sublattices on 
the pyrochlore lattice. 
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A, 


(0,1,0) 


(0,|,|) 


(- 1 -) 

V 2 ' ' 2 / 


(1,1,0) 


Bi 


(0,0,0) 


(o,|,|) 


(io,|) 


(1,1,0) 


bi 




(0,1,1) 


(1,0,1) 


(1,1,0) 



TABLE II: Primitive vectors (bi) and four basis vectors on 
each A-site (Ai) and B-site (Bi). 



Kondo-like coupling between A and B-sites 



In this section, we derive the explicit form of the 
Kondo-like coupling between the A-site non-Kramers 
doublet and the B-site Kramers doublet. As we argued 
in the main text, the Ising component of the A-site pseu- 
dospin Tj couples to their six neighboring B-site mag- 



netic moments ji 



2 L *P' 



whereas, the A-site pla- 



nar components rj couple to the B-site electron density 

rii = c-Cia. Hence, we separate the former and the latter 
cases and derive the Kondo-like coupling term allowed by 
lattice symmetry. First of all, let's consider the former 
case where t? couples to if. We note that this magnetic 
coupling term is already derived in Ref l28l 



TTZZ 

H AB 



E c i [til M + 3i (r~h + 63) + $ (r + h- b 2 ) + ft (r - b 2 + 63) + ]l(r) + f 3 (r + h- b 2 )} T *(r) 



+ {io ( r + b i) + Jo (r + h-b 2 + b 3 ) - f 2 - (r + h- 63) - j* (r + b 1 -b 2 + 63) 

- jl (r + 61) - if (r + h- 6 2 )K(r) + {f (r + h) + f (r + 63) 

- 3'(r) ~ fi (r + h) - f 3 (r) - Jf (r + MKM + {f Q (r + b 3 ) + f (r + b 1 -b 2 + b 3 ) 

- jl (r + b 3 ) - jf (r-b 2 + b 3 ) - f 2 (r-b 2 + b 3 ) - f 2 (r + b l -b 2 + 6 3 )}r|(r) 

+ c 2 [{jf (r) + jf (r-b 2 + b 3 ) + jf(r) + 3I (r - b 2 + 63) + jl (r + h- b 2 ) + if (r - b 2 + 63) 
+ 3 2 (r + h- b 2 ) + jl (r-b 2 + 63) + if (r) + if {r + h- b 2 ) + f 3 (r) + 3% (r + h~ b 2 )}r*(r) 
+ {-J y a (r + bi) - jl (r + h-b 2 + b 3 ) - j* (r + h) - j 2 (r + h - b 2 + b 3 ) 
+ jl (r + h- b 2 ) + jl (r + h-b 2 + 63) - jf (r + h- b 2 ) - j\ (r + h - b 2 + b 3 ) + if (r + h) 
+ if (r + h- b 2 ) - jl (r + h) - if (r + h- b 2 )}r*(r) + {-j* (r + h) - if (r + b 3 ) 

- 30 (r + h) - jl (r + b 3 ) - jf (r) - jf (r + b 3 ) + jf(r) + jf (r + b 3 ) 

+ jUr) + Jl (r + h) if (r) - if (r + h)}rl(r) + (r + 63) - if (r + h-b 2 + b 3 ) 

- io (r + b 3 ) ~ jl (r + h-b 2 + 63) - jf (r + b 3 ) - jf (r - b 2 + b 3 ) + jf (r + 63) + jf (r - b 2 + 

- if (r-b 2 + 63) - jf (r + h-b 2 + 63) + if (r-b 2 + 63) + if (r + 6 X - b 2 + b 3 )}ri(r) 



63) 



(6) 



Here, r?(r) is the Ising component of the pseudospin at 
sublattice i (— 0, 1,2,3) in unit cell labeled by r which 



is pointing along local Zi direction. (See Table [T] for their 



7 



local coordinate frames for the four sublattices on A-site 
pyrochlore lattice) In a similar way, j^{r) is y, (= x, y, z) 
component of B-site magnetic moments in global cubic 
coordinates. Primitive vectors (b{) and four basis vectors 
on each A-site (Ai) and B-site (Bi) are shown in Table [TT| 
There are two independent parameters, c\ and c 2 which 
can not be determined by symmetry grounds. However, 
two independent parameters can be easily understood as 
follows. For every given site i (on the A-site), there are 
one local Zi axis and its basal plane formed by x^ and 
j/i.(See Table [i]) The magnetic coupling can be written as 
T i(^jjj + c ^7+ c jjj) assuming that jj is aligned along 
their local axis. For its basal plane, we have freedom 



to choose Xj and yj, which leaves cj = c^. Once the 
coupling between A-site i and one of their six neighboring 
B-site j is determined, the coupling with all the other five 
neighbors are determined by lattice c 3 rotation and ah 
mirror symmetries, with two independent parameters Cj 
and = <?y 

On the other hand, the planar components of the A- 
site pseudospin are time reversal invariant and they do 
not couple to B-site magnetic moments but couple to the 
B-site electron density ni = c ia Ci a . In this case, there is 
only one independent parameter and one can write the 
coupling term as following 



h ab - c 3 



E 



[ r o~ ( r ){ uj2 ' n 2 (h -b 2 +r) + uj 2 n 2 (-62 + b 3 + r) + ujn 3 (61 - b 2 + r) + ujn 3 (r) 



+ «i (-62 + b 3 + r) + rii(r)} + T+(r){uj 2 n 3 (61 + r) + uj 2 n 3 (61 - b 2 + r) 
+ ujn 2 (61 - b 2 + r) + ljji 2 (bi - b 2 + b 3 + r) + n (61 + r) + n Q (bi - b 2 + b 3 + r)} 
+ T 2 ( r ){ uj2 ' n o {bi + r) + uj 2 n a (63 + r) + uni (b 3 + r) + n 3 (61 + r) + wni(r) + n 3 (r)} 

+ T 3 ( r ){^ 2 "i (^3 + r) + 0J 2 ni (—62 + 63 + r) + ujn Q (b 3 + r) + uono (61 - b 2 + b 3 + r) 

+ n 2 (-62 + b 3 + r)+n 2 {b x - b 2 + b 3 + r)} + h.c 



(7) 



Extra phases l,w = e l27r ^,uj 2 are present depending on 
the different sublattices due to their local coordinates 
on the A-site pseudospin. Combining H^ B and H^ B , 
the Kondo-like coupling with three independent coupling 
parameters c\,c 2 , c 3 is written as 



In the presence of magnetic ordering on the A-site, we 
take T s (r) (r s (r)) which acts as an effective magnetic 
field or chemical potential at every B-site itinerant elec- 
trons. 



H AB {c!,c 2 , c 3 ) = Hf B (ci, c 2 ) + H^ B (c 3 ). (8) 



